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1. Introduction 



1.1. Orientation 

Fusion rules yield the number of independent couplings between three given primary 
fields in conformal field theories. We are interested in fusion rules in unitary conformal field 
theories that have a Lie group symmetry, that is, those whose generating spectrum algebra 
is an affine Lie algebra at integer level. These are the Wess-Zumino-Witten models [1,2]. 
Primary fields in these cases are in 1-1 correspondence with the integrable representations 
of the appropriate affine Lie algebra at level k. Denote this set by Pj'^^ and a primary field 
by the corresponding affine weight A. Fusion coefficients " are defined by the product 



In the infinite-level limit and for fields with finite conformal dimensions (i.e. finite 
values of the finite (non-affine) Dynkin labels - see below), the purely affine condition on 
weight integrability is relaxed and the primary fields are solely characterized by their finite 
part, required to be an integrable weight of the corresponding finite Lie algebra. Recall 
that a weight A is characterized by its expansion coefficients in terms of the fundamental 
weights LOi 

r 

A = J]Aia;i = (Ai,...,A^) (1.2) 

i=l 

where r is the rank of the algebra. The numbers Aj's are the Dynkin labels. The set of 
weights with nonnegative Dynkin labels (the integrable weights) is denoted by P+. 

Similarly, in the infinite- level limit, the fusion coefficients reduce to tensor- product 
coefficients: 

hm A/-{')^=A6./. (1.3) 

fe— »oo '^A' 

where H\^^ is defined by 

\®^l=Y. -^a/^^ (1.4) 

By abuse of notation, we use the same symbol for the highest weight and the highest- weight 
representation. Notice that 

a6.m"=-A6.m-* (1-5) 
1 



where f* denotes the highest weight of the representation conjugate to that of f. Equiv- 
alently, Afx^u* gives the multiphcity of the scalar representation in the triple product 
X <S> n <^ i^* ■ 

This paper is devoted to a detailed analysis of tensor products from the point of view 
of generating-function methods, preparing the ground for an analysis of fusion rules along 
the same lines. This will be the subject of a following article. 



1.2. Tensor-product generating functions 

A tensor-product generating function codes the information for all the tensor products 
of a given algebra in a single function. This function is first defined in terms of an infinite 
series as follows. We introduce the dummy variables Lj,Mj,A^j with i = l,---,r and 
multiply the tensor product coefficient N'xfj," by L^M^N^ where 

L^ = L^---L^- (1.6) 

and sum over all integrable values of A, v: 

G{L,M,N)^ J2 Mxf.^'L^Mi'N'' (1.7) 

G can generally be expressed in terms of a rather simple closed function of its variables. 
For instance, for stt(2), it reads 

^' ^) = (l-LM)il-LN)il~MN) 

The above generating function contains all the global information concerning the 
tensor products in su{2). An example of global information that can be extracted from 
such an expression is the number of couplings having a particular property. Suppose that 
we are interested in the total number of triple products with \i + ni + Pi = 2n for a given 
positive integer value of n. We can rescale the three dummy variables by a factor s, giving 

G{L,M,N;s)= ^ 



(1 - s^LM){l - s^LN){l - s^MN) 
oo (1.9) 

n=0 Ai,^ii,i/i>0 



We then set L = M = N = 1 and read the desired information from the coefficient of s 
in the expansion: 

^ oo 
^ ' n=0 Ai+/[ii+z/i=2n 

that is 

Ai+/ui+i/i=2n 

Furthermore, the generating function contains local information. A simple example 
of useful local information is the value of an explicit tensor-product coefficient. Given a 
closed expression for the generating function, explicit tensor-product coefficients can be 
read from its Taylor series expansion. 

This in turn shows that another example of basic global information that can be 
deduced from a generating function is the integrality as well as the positivity of the tensor- 
product coefficients. In fact, tensor- product generating functions G can always be written 
in a manifestly positive form. 

More importantly, from our point of view, is that in the context of fusion rules, the 
construction of the simplest generating functions led to the discovery of the notion of 
threshold levels [3]. Moreover, as shown in the sequel paper, setting up a fusion generating 
function is a way to obtain explicit expressions for these threshold levels. Our new approach 
to fusion-rule generating functions, which originates from the generalization of techniques 
developed in the present paper on tensor products, leads to a further new concept, that of 
fusion basis. 



1.3. Overview of the paper 

The present article is organized as follows. We start by explaining in detail the con- 
struction of tensor-product generating functions for finite Lie algebras. The first construc- 
tion which is presented is the character method developed by Sharp and his collaborators 
(section 2). Here the starting point is the construction of the generating function for the 
characters of all the irreducible representations which serves as the building block for the 
construction of the tensor-product generating functions. Simple examples of generating 
functions are worked out in details (stt(2), sm(3), sp{A) and osp(l,2)). We also introduce 
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and exemplify the concept of a model for a generating function. This is an algebra whose 
Poincare series is the generating function under study. Models allow us to make use of 
powerful techniques from computational algebra to find generating functions as described 
below. 

Although it is conceptually very simple, the character method is strongly limited by 
its inherent computational difficulties: the disproportion between the simplicity of the 
resulting form of the generating function and the intermediate calculations is enormous. 
This motivates our alternative approach to the construction of tensor-product generating 
function. It is based on the reformulation of the problem of calculating tensor products in 
terms of the solution of a set of linear and homogeneous Diophantine equations (cf. section 
3). The Diophantine-equation reformulation of the problem is equivalent to its expression 
in terms of counting the number of integer solutions in a convex polytope and the polytope 
equations appropriate to any classical Lie algebra have been described by Berenstein and 
Zelevinsky. For su{N), these inequalities are equivalent to the Littlewood- Richardson al- 
gorithm, presented in section 4. The elementary solutions of these Diophantine equations 
represent "elementary couplings" . For sp(4) , the use of the Berenstein-Zelevinsky inequal- 
ities to obtain the elementary couplings and their relations (cf. the analysis of section 6) 
is new. 

The essential difficulty in constructing a generating function, given a set of inequalities 
is not finding the elementary solutions. For examples of the size we consider efficient 
algorithms exist to solve this problem. The key difficulty is related to the numerous 
relations - just alluded to - that exist in general between the elementary solutions. From 
the Diophantine-equation point of view, the decomposition of a solution may not be unique 
because different sums of elementary solutions could yield the same result. 

These redundancies must then be eliminated by forbidding one of the two sums of 
elementary solutions occurring in each relation. However, the procedure for doing this 
when there are more that one relation is somewhat tricky: we cannot select at random a 
sum to be forbidden from each relation and be sure that no further relations will arise. 
This is the problem of finding a compatible set of forbidden couplings. To solve this 
problem we first "exponentiate" it: given a solution a = (ai, . . . ,ak) to our system of 
linear Diophantine equations, we introduce formal variables Xi, . . . ,Xk and consider the 
monomial ^ . . . X^'^ . The linear span, R, of all such monomials is a model for the 
generating function for the solutions to the original set of linear Diophantine equations 
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(see section 5). To find the Poincare series of R we use the theory of Grobner bases. 
Given a set of elementary couphngs, or in this case generating monomials, Grobner basis 
techniques provide a method for generating a complete set of relations. In addition, and 
more importantly, they provide an algorithm for specifying a compatible set of "forbidden 
couplings". In other words one can calculate the Poincare series of R and hence the 
required generating function, directly from the set of relations provided by the Grobner 
basis method. It is then shown how this machinery applies to the construction of tensor- 
product generating functions. This is the subject of section 5. 

A trivial modification of the description of tensor products in terms of inequalities 
follows by introducing extra positive labels that transforms the inequalities into equalities. 
For su{N), this gives rise to a remarkable graphical construction, the famous Berenstein- 
Zelevinsky triangles. These are introduced in section 7. We also discuss the analogous 
construction for sp{4), whose diagrammatic representation is new. But the main interest 
of these reformulations is that it yields a simple and systematic way of obtaining the 
elementary couplings from the construction of a vector basis. A projection operation 
yields the elementary couplings. In the process, we get a new way of constructing the 
generating functions. This vector-basis approach to the construction of tensor-product 
generating functions is illustrated with su{3) and sp(4). 

1.4- What is new? 

In part, this article is a comprehensive review of existing techniques for constructing 
tensor-product generating functions. In addition, the efficiency and the power of the 
different methods of construction are compared. We also, from a practical point of view, 
discuss the connections between commutative algebra and the computation and properties 
of generating functions. On the mathematical side, these links have been developed by 
Stanley [4]. In this article, emphasis is placed on the concept of a model of a generating 
function, which is an algebra whose Poincare series is equal to the given generating function. 
This allows us to make use of powerful techniques in computational commutative algebra 
including Grobner bases, which were first used in this context in [5]. 

In this way we arrive at new derivations of previously known results. But those 
pertaining to sp(4), and in particular the use of the Berenstein-Zelevinsky inequalities as 
a way to fix unambiguously the linear relations between the "elementary couplings", are 
new and should be stressed. However, the main achievement is a sound reformulation of 
the problem that is well-adapted to its extension to the fusion case. 
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2. Generating-function for tensor products: the character method 



2.1. The character method for the construction of the tensor-product generating function: 
the su{2) case 

The method developed by Sharp and coUaborators for constructing generating func- 
tions for tensor products is based on manipulations of the character generating functions 
[6]. Although simple in principle, these manipulations become rather cumbersome as the 
rank of the algebra is increased. To illustrate the method, we will work in complete detail 
the simplest example, the su{2) case. 

The first step is the derivation of the character generating function. The Weyl char- 
acter formula for a general algebra of rank r and a highest-weight representation A is 

XA = ^ (2.1) 

where p is the finite Weyl vector, p = Yli=i '^ii ^^'^ where the characteristic function ^ is 
defined as 

Ca+p = (2.2) 

wew 

where e{w) is the signature of the Weyl refiection w and W stands for the Weyl group. 
For stt(2), the Weyl group contains two elements: 1, si. With 

X = e'^i (2.3) 

the su{2) characteristic function ^ for the representation of highest weight muji = (m) is 

U+i = x^+' - X— 1 (2.4) 

The character reads then 



Xm = ^^T— = , -2 = + + --- + x-^ (2.5) 

tAJ iAj -L tXj 

The character generating function xl is obtained by multiplying the above expression by 
L™ where L is a dummy variable, and summing over all positive values of m: 

oo oo 

Xl{x) = L^xm = n E L'-i^'^^' - ^"™"') 

X — X ^—^ 

(2.6) 

1/1 x-^ \ 1 



l-x-'^\\-Lx l-Lx-^J {I - Lx){l - Lx-^) 
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By construction, the character of the highest weight (m) can be recovered from the power 
expansion of xl as the coefficient of the term L"^. The characteristic generating function 
is defined by 

XL{x) = ^f (2.7) 



and it reads 



X — X ^ X X ^ 



^^(^) - (1 _ Lx){l - Lx-i) ~ \ -Lx 1 - Lx-^ ^^-^^ 
the last form being the one that results directly from (2.5). 

The tensor product of two highest-weight representations can be obtained from the 
product of the corresponding characters: 

XmXn = ^N'mJ Xl (2-9) 

I 

This information can be extracted from the product of the corresponding generating func- 
tions. We are thus led to consider the product Xl(3^)xm(3^)- To simplify the analysis of the 
resulting expression, notice that the information concerning the representations occurring 
in the tensor product is coded in the leading term of the character, i.e., the term x"^"*"^. 
The rest of the representation is easily reconstructed by the application of the Weyl group 
and the action of the ladder operator. Actually, to insure that every positive power of x 
singles out a highest- weight representation, we can multiply both sides by To read off 
these terms, we can focus on the terms with strictly positive powers of x in the product 
XLix)xM{x)^Q{x)- If we want to read off directly the Dynkin label of the representations 
(and not their shifted value), it is more convenient to divide by x before doing the projec- 
tion, now restricted to the nonnegative powers of x. The truncation of an expression by 
its negative powers of x will be denoted by the MacMahon symbol [7] defined by 



oo 

X 



^ Cnx" = ^ c„a;" (2.10) 

~ — oo n>0 

When there is no ambiguity concerning the variable in terms of which the projection is 
defined, it is omitted from the symbol. 

We are thus interested in the projection of the following expression 

XL{x)xMix)iQ{x)x~^ = Xl{x)^m{x)x~'^ 



x-^ 



(1 - Lx)(l - Lx-^) \1-Mx 1 - Mx-^ 
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(2.11) 



For these manipulations, we use systematically the following simple identities: 

1 _ 1 / 1 Bx-'^ 

(1 - Ax)(l - Bx-^) ~ (1 - AB) \1-Ax ^ 1 - Bx-^ 

^ ( Ax 1 , .„-.„x 

+ ^-rr (2-12) 



{1- AB) Vl - via; 1 - Bx 



(1 - AB) \\- Ax^ \- Bx-^ ^ 

Which one is used is dictated by the context (a good choice often simplifies the subsequent 
analysis, the last one being somewhat more algorithmic). 

We thus have two terms to analyze. The first one is 

' 1 (2.13) 



(1 - Lx){\ - La;-i)(l - Mx) (1 - La;)(l - LM) \\-Mx 1 - Lx'^ 
The first part is not affected by the projection. The second one can be written as 

Lx~^ Lx~^ ( Lx 1 \ /oi/i\ 

(1 - Lx){l - LM){1 - Lx-^) " (1 - LM)(1 - L^) \T^^ ^ 1 - Lx'^ J ^ 

The second part of the above expression contains only negative powers of x and can thus 
be ignored. The first part is unaffected by the projection. We have thus, for the first term 
of (2.11) 

1 1 / 1 \ 

^{1-Lx){l-Lx-^){1-Mx)~ {1-Lx){l-LM)\1-Mx'^ l-L^J ^' ^ 

The other term that needs to be considered is 

_2 

X 

n 



> (1 - Lx){l - Lx-^){1 - Mx-^) 

x"^ / 1 Mx 

7 TTT^ ^TTT Z ^ + 



> {1- Lx-'^){1- LM) \1- Lx l-Mx-i 
x~^ 

n 



> (1 - Lx-i)(l - LM){1 - Lx) 



> (1-LM)(1-L2) {l-Lx^ l-Lx-i J ^^""^^^ 

= n 



> (1 - LM)(1 - L'^){1 - Lx) 
Lx-^ ( 1 



(1-LM)(1-L2) \\-Lx 

L2 

{\ - LM){\ - L'^){\ - Lx) 



- 1 



Subtracting (2.16) from (2.15), we find that 



Replacing x by A^, we thus get 

G'^'C^hL.M.N) = (2.18) 

^ ' ' ^ {1- LM){1- LN){1- MN) ^ ' 

as the generating function for su{2) tensor products. To read off the representations (n) 
that occur in the tensor product {£) ® (m), we expand the generating function (as a Taylor 
series) and collect all terms A^" that are multiplied by L^M^. All such values of n are the 
representations we are looking for. This is simply a restatement of the following relation 
between the tensor- pro duct coefficients MurrT" and the generating function: 

G{L,M,N)= J2 J^inJ'L^M'^N'' (2.19) 



2.2. The abstract setting: Poincare series, elementary couplings and relations; defining a 
model 

As we shall see it is frequently useful have a model, R, for a generating function 
G{Xi, . . . , Xk) such as (2.18). By this we mean a commutative Q-algebra with an identity, 
graded by N'^, 

-R = ©QgNfe-Ra 5 RaRp ^ Ra+I3 (2.20) 

and such that its Poincare series i 



F{R) = dimQ(i?«)X^ 



satisfies 

F{R) = G{X^,...,Xk). (2.21) 

For example, for (2.18), with Xi = L, X2 = M, Xs = N, we can take R = Q[Ei, E2, Es] 
(in fact all our examples R is either a subring or quotient of a polynomial) with the 



^ Such series are also called Hilbert series. 
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grading of -Ei, -E2 and E3 being (1, 1, 0), (1, 0, 1) and (0, 1, 1). The homogeneous subspaces 
are spanned by i^f £^2^3, a, 6, c e N with grade {a + b, a + c,b + c) and so 

F{R)= J2 X^+''X^+'=X^+'' = G'''^^\X^,X2,X3) (2.22) 

(a,6,c)eN3 

as required. 

If R is generated by elements Ei, . . . , Eg and is a model for a generating function 
G for tensor products (or fusion products) then we call Ei, . . . , Eg a set of "elementary 
couplings" for G. 

It should perhaps be stressed that a priori the variables Xi, . . . , Xk and Ei, . . . ,Eg 
are unrelated. We shall refer to the £"s as model variables and the X's as grading variables. 
If the grading vector of £'j is a* , z = 1, . . . , s then there is an associated monomial in the 
grading variables: X°'\ for which we will use the notation g{Ei). For example in the above 
example we have g{Ei) — X\XlX^ = LM. However, to avoid tedious repetition when 
writing down generating functions we shall often write, for example, 1/(1 — £^1) rather than 
1/(1 — fif(£'i)). In all such cases where model variables appear in a generating function they 
should be replaced by the corresponding monomial in the grading variables. The utility 
of this convention can be seen by examining some of the examples of generating functions 
contained in the following sections and comparing them with the much more unwieldy 
expressions obtained by everywhere replacing the model variables by grading variables. 

In the case of tensor products we use the notation : g{E) : product" to denote 
a set of elementary couplings with their "exponentiated" grading and the corresponding 
term in the tensor product. So in the example above we would write: 

Ei-.LM: (1)®(1)D(0), 

E2:LN: (1) ® (0) D (1), (2.23) 
E3 : MN : (0) (1) D (1) 

Having made the distinction between grading and model variables, it should be noted 
that there are cases where we can identify the model as a ring generated by monomials in 
the grading variables. So in the above example we could define Ex = LM, E2 = LN and 
Es = MN and take the model for our generating function to be the subring of Q[L, M, N] 
generated by Ei, E2 and £'3. However, it is not always desirable, or even possible, to make 
this identification. 
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As we shall see, identifying a generating function with the Poincare series of a model 
will allow us to link the generating function more explicitly with combinatorial rules for 
calculating tensor product (and other) multiplicities. It will also allow us to make use of 
powerful algebraic techniques for computing Poincare series. 

In the example above we are given a generating function from which it is easy to find 
a model. Of course in general we will start with a model from which we construct the 
required generating function. We close this section with two examples of how models for 
the su{2) character generating function can be constructed. 

The first method, which has been exploited by Sharp et al (see [6]) to construct 
character generating functions, amounts to finding an algebra R which is a module for the 
Lie algebra su{2) and such that as an su{2) module R is isomorphic to (Bi>iVi where Vi is 
the irreducible su{2) module of dimension i. 

In this case we can take R = Q[p, q] with the generators of su{2) being given by 
differential operators: 

d d d d 

h=p-^ -q-^, x- = q—, x+=p— (2.24) 
op oq op oq 

The su{2) highest-weight vectors are p*, z > and a basis of the irreducible submodule 
of dimension i is just given by the monomials of degree i m. p and q. We can give R an 
grading by taking the degree of p to be (1, 1, 0) and of q to be (1, 0, 1). Here the first 
grading index specifies the representation while the other two refer to a particular weight. 
As i? = Q[p, q] the Poincare function for R is, 

1 



(l-p)(l-g) 

with the understanding, as explained above, that p and q should be replaced by the corre- 
sponding expression in terms of the grading variables. Let us denote these grading variables 
here by L (which exponentiates the representation index) and x, y (exponentially related 
to the weights) . The Poincare function reads then 

(2.25) 



{1-Lx){l-Ly) 



Of course in this case the distinction involved and the usefulness of our convention is not 
immediately apparent. 
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Note that in this example we could take a 1? grading with the grades of p and q being 
(1,1) and (1,-1). However in general allowing negative weights means we have to be much 
more careful when taking Laurent series expansions and so we avoid this complication when 
possible. This grading is related to a single variable description of the weights • • • , a;~*. 
Setting y — gives the corresponding form of the generating function. 

Another way of constructing a model for the weight generating function, which makes 
more natural the grading, is to observe that the complete set SU{2) weight vectors 
of finite dimensional irreducible su{2) modules are in 1-1 correspondence with one-rowed 
Young tableaux. If the Young tableau has c boxes filled with a I's and b 2's then there is 
a constraint 

a + b-c=0, a,6,c>0 (2.26) 

and so the solutions to this linear Diophantine equation are in 1-1 correspondence with 
the complete set of SU (2) weight vectors. Thus to find a model for the weight generating 
function it is sufficient to find a model for the solutions to (2.26). In this example it is 
not difficult to see that every solution to this equation is a linear combination (with non- 
negative coefficients) of the two fundamental solutions: (a, 6, c) = (1,0,1) and (a, 6, c) = 
(0, 1, 1). Let R be the subring of Q[^, B, C] generated by the monomials Ei = AC, E2 = 
BC. Considering the exponents of the monomials Ei and £'2, we see that the monomials in 
R correspond to the solutions of (2.26) and hence taking the natural grading on R ensures 
that the Poincare series of R is the generating function for the solutions to (2.26) and 
hence is the required generating function. In this example there are no relations between 
El and E2 and so R is isomorphic to the polynomial ring in two variables (as expected) 
and so the Poincare function is once again (with A — > x, i? — > y, C — > L): 

{l-Lx){l-Ly) ^^-^^^ 

In all these simple examples we have found a model which is a polynomial ring. In 
general this will not be the case. This is illustrated in the next section. 

2.3. Multiple su{2) tensor products 

In order to illustrate in a rather simple context the occurrence of relations between 
elementary couplings, we will consider a slight generalisation of the previous problem 
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namely, finding the multiplicity of a given representation ( in the triple product A ® /U ® 
u. In terms of character generating functions, this amounts to considering the product 
Xl{x)xm{x)xn{x) D Xp{x), or equivalently, Xl{x)xm{x)^n{x)x~'^ D ^p{x)x~'^. The left 
side is then projected onto positive powers of x. We are thus led to consider 

1 x-^ \ 

> (1 - Lx){l - Lx-i)(l - Mx){l - Mx-i) Vl - A^x ~ 1 - Nx'^ ) 

The projection of each term is worked out as previously and the resulting expression is 
found to be, with x replaced by P: 

G{L, M, N, P) = (1 _ _ MP)(1 - NP){1 - LM)(1 - LN){1 - MN) ^^'^^^ 

This is the sought for generating function. Here we would like to have a model with 6 
elementary couplings corresponding to the terms in the denominator of the generating 
function: 



El : 


: LM : 


(1)®(1)®(0) D (0) 


i?2 : 


: LN : 


(1)®(0)®(1) D (0) 


Es : 


: LP : 


(1)®(0)®(0) D (1) 




: MN : 


(0)®(1)®(1) D (0) 




: MP : 


(0)®(1)®(0) D (1) 


Eq : 


: NP : 


(0)®(0)®(1) D (1) 



(2.30) 



and there must be a linear relation between the following products (signalled by a term in 
the numerator) which have grading LMNP: @ 

EiEq, E2E5, EsE4 (2-31) 



It is not difficult to see that a model is given by Q[ei, 62, 63, 64, 65, ee]// where Ei = 
Ci + /, z = 1, . . . , 6 and / = (aeiCe +66265 + C6364) is the ideal generated by the polynomial 
a6i66 + 66265 + C6364 with a, 6, c G Q not all zero. In particular the model is far from 
unique. As we shall see later, particular methods of construction will select one particular 
model. 

Relations between products of elementary couplings are often called syzygies in the physics 
literature (see in particular [6] and related works). However, the proper mathematical meaning 
of a syzygy is somewhat different and for this reason we stick to the more correct terminology of 
relation. 
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Before leaving this example, we would like to rework it from a different point of view, 
as an illustration of the 'composition' technique of generating functions. Let G{L, M, R) 
describe the tensor product corresponding to XlXm ^ Xr and similarly let G{Q, N, P) 
correspond to XQXn 3 Xp- We are interested the product Xl{x)xm{x)xn{x) D Xp(a;) 
but treated from the product of the two generating functions G. We thus want to enforce 
the constraint R = Q in the product G{L, M, R)G{Q, P). The ideai is to multiply this 
product by (1 — R~^)~^ and, in the expansion in powers of R and Q, keep only terms 
of order zero in both variables: with an obvious notation we have 

h?lG{L,M,R)G{Q,N,P) ^ 



1-Q-^R-^ 

= §§ 5^ M) R^ Yl ^rn{N, P) Q"^ (2.32) 

n m a 

= 5^Ap(L,M)Sp(L,M) 

V 

which is manifestly equivalent to considering 

VlG[L,M,x)G{x-^,N,P) (2.33) 

With the explicit expressions for the generating functions, we have thus 
X \ 1 

n I t (2 341 

= (1-La;)(l-Mx)(l-LM) (l-Pa;-i)(l-A^a;-i)(l-A^P) ^' ' 

A brief and by now standard analysis yields directly the generating function (2.29). 

2.4. The osp{l,2) case 

The simplest example after su{2) is that of the superalgebra osp(l,2). Very little 
information is needed about superalgebras for the study of the osp(l, 2) representations. 
We only need the fact that the highest weight osp(l,2) representations {m} decompose 
into a direct sum of two su{2) representations; (m) © (m — 1/2) [8]. The character is thus 

XW = — -^-^ (2-35) 



^ This trick is used in different references in [6], mainly in relation with the construction of 
generating functions for branching functions. 
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It is a simple exercise to check that the generating character function is 
The tensor-product generating function is found to be [9] 



\2 



^ ' ' ^ {l-LM){l-LN){l-MN){l-LMN) 

1 + LMN 
~ (1 - LM)(1 - LN){1 - MN) 

An underlying model would thus have four elementary couplings: 

E^-.LM: {1}®{1}d{0} 
E^-LN: {1}«){0}d{1} 
^3 : MN : {0} ® {1} D {1} 
^4 : LMN : {1} ® {1} D {1} 



(2.37) 



(2.38) 



and the numerator indicates a linear relation between products of elementary couplings of 
degree (LMN)^ namely: 

E1E2E3 and Ej (2.39) 

A model is given by Q[ei, 62, 63, 64]// where Ei = Ci+I, i — 1, ... ,4 and I = (0616263 + 
664) is the ideal generated by the polynomial 0616263 + 664 with a, 6 G Q, but not both 
zero. 

2.5. The su{3) case 

The next example in complexity is su{3). With 

Xi = e'^' i = l,2 (2.40) 
the characteristic function for a representation of highest weight (m, n) = muji + nu!2 is 

t, , _ ^rn+l n+l _ -m-1 m+n+2 _ n+m+2 -n-1 
s(m,n) — -^-l •*'2 •*'2 "^1 "^2 

I ™n+l -m-n-2 I -m-n-2 m+1 _ „-n-l -m-1 
15 



The characteristic generating function is obtained by multiplying this result by U^L^ and 
summing over all positives values of m, n. The result is 

{1 - LiXi)(l - LiX^^X2){l - LiX2^){l - L2X2){1 - L2XiX2^){l - L2X^^) 

(2.42) 

The construction of the tensor-product generating function proceeds exactly as for 
su{2), but here it is much more complicated from the simple fact that there are two 
variables. The result is [6] 

G'^"(3)(Li, L2, Ml, M2, A^i, N2) = (1 - LiLsMiMsA^iA^s) 
X [(1 - LiN2){l - LiM2)(l - L2Mi)(l - LaA^i)]-' (2.43) 
X [(1 - M2iVi)(l - MiAr2)(l - LiMiiVi)(l - L2M2iV2)]-' 

From the denominator we see that it is natural to seek a model with eight elementary 
couplings: 

El : L1M2, E2 : LiiV2, Eg : MiN2 

E4 : L2M1, E5 : L2N1, Eg : M2N1 (2.44) 

Er : LiMiNi, Eg : L2M2N2 
The numerator indicates a relation between these elementary couplings that need to be 
taken into account to avoid over-counting multiplicities: there must be a linear relation 
between the following three products: 

E1E3E5, E2E^Eq and E^E^ (2.45) 

which are the three terms with grading L1L2M1M2N1N2. This is the only relation required 
and a model for this generating function is given by i? = Q[ei, . . . , eg]//, with Ei = 
ei + I,i = 1, . . . ,8 and / = (aeiCgCs -|- 6626466 + 66763) with 0616365 -|- 6626466 + 66763 7^ 0. 

The elements of R have the form m + I with m e Q[ei, . . . , eg]. However there is no 
canonical way of choosing the representatives m. Take for example the case a = 6 = c = 1. 
(Usually we will construct a model for our generating function as explained above and this 
construction will fix the values of a, b and c). In R we have EiE^E^ = —{E2E4EQ + EjEg) 
and so we can take as a basis for R the set of (equivalences classes of ) monomials which do 
not contain the product E1E3E5. In this case we say that we have chosen to make E-^E^E^ 
a 'forbidden product'. Similarly we can forbid the products E2E4EQ or E-jEg. As we shall 
see later, the choice of forbidden products corresponds to a choice of term ordering. 
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We can write the generating function in three different ways, each making manifest the 
fact that one of the above products never appears. Note once again we use the convention 
that Ei should be replaced by g{Ei) in these expressions: 



G 



\i^l,3,5 / ^ 



El 



(l-EOa-Es) 



E3 _|_ E5 



{1-Es){l-Ei) {1-E5){l-Es) 

E' 




(2.46) 



{1-E^){1-E^) {1-Ee){l-E2) 



It is clear that in expanding the first form, we will never encounter a term corresponding 
to a product of the three factors EiE^E^. Similarly no product E2E4EQ corresponds to 
a term in the second form, while the last expression amounts to forbidding all factors 
containing E^Eg. Therefore, although G is unique, its expression in terms of the Ei^s is 
not, because the forbidden couplings may be chosen in different ways. 

Before leaving this example, let us mention another way of constructing the generating 
function using the idea of 'composition' of generating functions described previously. This 
uses the Giambelli formula that expresses a general representation in terms of a difference 
of products of representations with a single non-zero Dynkin label, i.e., 

(/Xi, //2) = (//I + ^2, 0) ® (/i2, 0) - (//I + //2 + 1, 0) ® (/i2 - 1, 0) (2.47) 

First consider the generating function: 

G{Li,L2,Mi,Ni,N2) (2.48) 

which is the generating function for products of the form: (Ai,A2) (//i,0). From this 
generating function we form: 

H{Li,L2,Mi,M2,Ri,R2) = ?lG{Li,L2,Mi,Ni,N2)G{N^\N^\M2,Ri,R2) (2.49) 
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which is the generating function for products of the form 

(Ai,A2)® (mi,0)® (//2,0) (2.50) 
Note that the generating function for products 

(Ai,A2)(8)(//i + l,0)®(//2-l,0) (2.51) 
is HM2Mi^ and so, by (2.47), the generating function for products (Ai, A2) (X> (a^i, ^^2) is: 

n[H- HM^M:^^] (2.52) 

The coefficient of Mf^M^*^ is the multiphcity of the representation with Dynkin labels 
(fii — //2) A*2) in the product 

(Ai, A2) ® [(/xi, 0) ® (/X2, 0) - (/xi + 1,0)® (/X2 - 1, 0)] (2.53) 

To change to variables which carry the Dynkin labels we make the substitution M2 1— > 
M2M]~"^, so that Ml now carries the first Dynkin label. This introduces negative powers 
of Ml, corresponding to products 

(Ai, A2) ® [(//I, 0) {112, 0) - (//I + 1, 0) {112 - 1, 0)] (2.54) 

with fii < fi2, which are not required. So we must keep only positive degree terms in Mi 
to obtain the final generating function. 

2.6. The sp{4) case 

As a final example, consider the sp(4) case. With the Xi defined as in (2.40), the 
characteristic function is found to be 



m+l n+l — m — 1 m+n+2 ™'^+'ti+5 — n — 1 , m+2n+3 — m— n— 2 



S(m,n) — "^1 -^2 -^l "^2 "^1 "^2 



I „-m— 2n-3 n+m+2 _ m+l -m— n-2 _ -m—2n-3 n , — m-1 -n— 1 
~l~ tJu 2 2 *^2 1*^2^ 2 

and the characteristic generating function is 

1 



(2.55) 



/ 



(1 - Lia;i)(l - Lia;ia;2 ^)(1 - L2X2 ^)(1 - -£'2iCi ^3^2) 
l + i^2 (l + L2)Lia;i 

(1 - L2X\X2^){1 - L2X2^) (1 - LiXi)(l - L2XIX2^) ^^"^^^ 
LiX^^X2 



(1 - Lia;i)(l - Lia^i ^0:2) 
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From this characteristic generating function, we construct the character generating func- 
tion and then we can proceed to the tensor-product generating function. This is again 
extremely cumbersome. The result is [10] 

G'P^^^iL,,L2,M^,M2,N^,N2) 

= [(1 - MiAri)(l - LiNi){l - LiMi)(l - M2iV2)(l - L2N2){1 - L2M2)]-' 
, 1 ^ L2M2NI 



(l-L2MiiVi)(l-L2M2Ar2) (l-L2MiAri)(l-L2M2iV2) ^^.57) 

LlM^NiN2 L1M2N1 



+ 



(1 - LiM2iVi)(l - LIM2N2) (1 - LiM27Vi)(l - L2M2NI) 
LIM2N2 L1M1N2 



+ 



(1 - LiMiiV2)(l - LIM2N2) (1 - LiMiiV2)(l - L2MIN2), 

From this expression, we read off the following list of elementary couplings (recall that the 
first variable is a model variable and then we write the corresponding monomial in the 
grading variables): 

Ai-.MiNi, A2:LiiVi, A3 : LiMi 

Si : M2iV2, B2 : L2N2, B3 : L2M2 

(2.58) 

Ci : L2M1N1, C2 : L1M2N1, C3 : L1M1N2 
Di : LIM2N2, D2 : L2MIN2, Ds : L2M2NI 

However, not all the products of the model variables can be linearly independent: there 
are linear relations between: 

QCj, AkDk, and AiAjBk 

DiDj, AlB.Bj, and BkCl (2.59) 

CiDi, AjBkCk, and A^BjCj 

for k a cyclic permutation of 1, 2, 3 (and repeated indices are not summed). (It is plain 
that the three sets of products found to be linearly related must have the same Dynkin 
labels.) A specific form of the generating functions amounts to a specific choice of a set of 
forbidden couplings among those that are related by a linear relation. 

2. 1. Technical remarks on the character method 

The character method provides a first principle approach to the construction of tensor- 
product generating functions. This is certainly its great virtue. However, the last two 
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examples indicate the essential limitation of this method for constructing generating func- 
tions for tensor products: the calculations are extremely complicated. Further progress 
requires the search for more powerful techniques. 

It should also be clear from the last example that in general the problem of finding 
a model and computing its Poincare series is not trivial even when we know the form of 
the generating function. Of course we shall be interested in the inverse procedure, that 
is, in constructing models in order to compute the associated generating function. This 
requires a systematic procedure for computing Poincare series and this will be provided 
by the theory of Grobner bases which we shall discuss shortly. 

3. Tensor-product descriptions 

3.1. The need for a tensor-product description 

As already mentioned, the fundamental limitation of the character method for con- 
structing tensor-product generating functions lies in the complication of the intermediate 
calculations associated to the projection to positive powers of the XiS variables of the 
different terms of the generating character products. The complication of these interme- 
diate steps should be contrasted with the relative simplicity of the resulting generating 
functions. This state of affair strongly suggests that there are much more efficient ways of 
obtaining these generating functions. 

It is clear that one major technical complication of the character method is that it 
starts at too fundamental a level, namely the character of the separate representations. As 
a result, we need to take care of the action of the Weyl group: this generates many terms 
and the projection of each term is rather complicated. But the fact that their final sum 
conspires to produce a rather simple result suggests that bypassing the use of the Weyl 
group would induce substantial simplifications.i 

One natural way to proceed is to start from a combinatorial description of the tensor- 
product rules. Such a description already takes into account the action of the Weyl group 
and encodes the various subtractions of the singular vectors. 

^ This continues to be the case when affine Weyl reflections are included for the fusion-rule 
generating functions, with the additional complication that extra projections are required to ensure 
that the representation labels are in the fundamental region of the affine Weyl group. 
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But how do we make the connection with the generating-function approach? The key 
is to find a combinatorial description which can be expressed as a set of hnear Diophantine 
inequahties. 

Given this set of inequahties, there is an algorithm, again due to MacMahonll, for con- 
structing a generating function. In this context, the generating-function method appears 
as a general approach to the solutions of a system of inequalities. This is particularly well 
illustrated in MacMahon's book [7]i. This provides then a direct route from the Diophan- 
tine inequalities to the generating function. This method is conceptually similar to the 
character method, except that the starting point is substantially closer to the end result. 
Let us illustrate MacMahon's approach with a simple example. 

3.2. MacMahon's theory for the solution of Diophantine inequalities 

We look for all the positive integer solutions of the inequahties (see e.g. [7] vol. 2 no 
356 p. 109) 

Xi > X2 Xi > Xs (3-1) 

To impose these constraints on a free series of the form 

J2^l'^2^^3' (3-2) 

we introduce extra parameters t and s as follows. To take into account the first inequality, 
we replace Xi by tXi and X2 by t~^X2 and project onto positive powers of t. Similarly, 
the second inequality is taken care by the replacements Xi sXi and X^ s~^Xs and 
projecting onto positive powers of s. This leads to 

t s 1 

> > il-tsX,){l-t-^X2){l-s-^Xs) ^^-^^ 

whose projections read 

1 - StXfX2X3 



(1 - tsXi){l - sXiX2)il - tXiXs){l - X1X2XS] 



(3.4) 



^ This is an adaptation of a method developed by Elhot [11] for the analysis of linear Dio- 
phantine equalities and for this reason the algorithm is often referred to as the Elliot-MacMahon 
method. 

^ See in particular vol. 2 section VIII. 
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We then set the auxihary variables s, t equal to 1 and obtain: 

1 — XfX2Xs 



(3.5) 



(1 - Xi)(l - XiX2){l - XiXs)il - X1X2XS) 

From this we read the four elementary solutions to the Diophantine inequalities under 
study: 



ai = (1,0,0), a2 = (1,1,0) 

as = (1,0,1), a4 = (1,1,1) 
with the ordering {xi,X2,xs) and the linear relation 



(3.6) 



a2 + as = ai + ot^ (3-7) 

3. 3. Seeking a road from elementary couplings to generating functions 

Although the description of tensor products via linear Diophantine equations is a more 
efficient route to finding the generating function than the character one, complications 
associated to the O projections remain a source of technical difficulty that severely limits 
the practical applicability of the method. 

A more powerful approach to our problem is to use the techniques of computational 
algebra. We start with a description of the tensor-product multiplicities as solutions to 
linear Diophantine inequalities. From these we find directly a model for the generating 
function.i 

In this approach, we thus work with a particular model of our generating function. 
But the advantage of this is that we can use Grobner-basis techniques, described below, 
to find the Poincare series of the model and hence the required generating function. 



4. The LR rule isu{N)) 

For su{N) tensor products, there is a particularly convenient description, which is that 
of the Littlewood-Richardson tableaux, supplemented by the stretched-product operation 
(defined below) [5]. 

^ This is roughly the inverse of the MacMahon's method which was originally conceived as a 
technique to generate the elementary couplings and their linear relations through the construction 
of the generating function. Here, the elementary couplings and their relations are first obtained 
and used as the input for the construction of the generating function. 
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Integrable weights in su{N) can be represented by tableaux: the weight (Ai, A2, • • • , Ajv-i) 

is associated to a left justified tableau of N — 1 rows with Ai + A2 H h Aat-i boxes in the 

first row, Ai + A2 + ■ • • + AAr_2 boxes in the second row, etc. Equivalcntly, the tableau has 
Ai columns of 1 box, A2 columns of 2 boxes, etc. The scalar representation has no boxes, 
or equivalently, any number of columns of N boxes. 

The Littlewood-Richardson rule is a simple combinatorial description of the tensor 
product of two su{N) representations A®//. The second tableau (//) is filled with numbers 
as follows: the first row with I's, the second row with 2's, etc. All the boxes with a 1 are 
then added to the first tableau according to following restrictions: 

1) the resulting tableau must be regular: the number of boxes in a given row must be 
smaller or equal to the number of boxes in the row just above; 

2) the resulting tableau must not contain two boxes marked by 1 in the same column. 

All the boxes marked by a 2 are the added to the resulting tableaux according to the above 
two rules (with 1 is replaced by 2) and the further restriction: 

3) in counting from right to left and top to bottom, the number of I's must always be 
greater or equal to the number of 2's. 

The process is repeated with the boxes marked by a 3, 4, • • •, with the additional rule that 
the number of z's must always be greater or equal to the number of z + I's when counted 
from right to left and top to bottom. The resulting Littlewood-Richardson (LR) tableaux 
are the Young tableaux of the irreducible representations occurring in the decomposition. 

These rules can be rephrased in an algebraic way as follows. Define to be the 
number of boxes i that appear in the LR tableau in the row j. The LR conditions read 



fc-i 



min{k,N 



1) 



^.-1 + E 



1 — 



E 



riii > 



i<A;<i<A^ j 7^1 



(4.1) 



and 



k 



k 




2<i<k< N et z<iV-l. 



(4.2) 
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The weight ji of the second tableau and the weight v of the resulting LR tableau are 
respectively given by 

N N-1 

^^ij^^N z = 1,2, ...,7V- 1 , 

(4.3) 

N-l min(j,iV-l) 

- + ^riij+i = Hij j = 1,2,...,N -1 . 

i=l i=l 

Hence, given three weights A, /U and z/, the number of positive integers solutions {riij} 
satisfying the above conditions gives the multiplicity N^^'^ of u in the tensor product 

A (8) 

The combined equations (4.1) and (4.2) constitute a set of linear and homogeneous 
inequalities. As described in [4] , the Hilbert basis theorem guarantees that every solution 
can be expanded in terms of the elementary solutions of these inequalities. 

As explained in section 3.2 we can construct a model for the solutions of the equations 
(4.1) and (4.2) by introducing new formal variables Ai, 1 < i < t where t is the total 
number of variables in (4.1) and (4.2). Then the subring of Q[Ai; 1 < i < t] generated 
by the monomials A" with a a solution of (4.1) and (4.2) provides the required model. 
This ring R will be generated by a finite set of monomials Ej 1 < j < s which we call 
elementary couplings corresponding to the elementary solutions of (4.1) and (4.2). Thus R 
is isomorphic to Q[ei, . . . , Cg]// under the mapping cj) : ei ^ Ei where / is an ideal. Each 
element of / corresponds, via the map 0, to a relation between the elementary couplings. 
As we shall see shortly, calculating a generating set of elements of / (or more particularly 
a Grobner basis) is the key step in our calculations. 

In the case of LR tableaux, there is a nice pictorial representation of the model R. 
Consider the set of formal linear combinations of LR tableaux with rational coefficients. It 
is given a ring structure by defining the stretched product of two LR tableaux (denoted by 
•) to be the tableau obtained by fusing the two tableaux and reorder the numbers in each 
row in increasing order. More algebraically, if we denote the void boxes of a LR tableau 
by a 0, so that 

AT-l 

noj=Y.\ j = l,2,...,iV-l (4.4) 

i=j 

we can characterize completely a tableau by the data {ny} with now z > 0. It is clear 
the set of numbers {n^} with z > 0, or equivalently, {Ai,ny} with z > 1, is a complete 
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set of variables for the description of the tensor products. Then, the tableau obtained by 
the stretched product of the tableaux {riij} and {n[j} is simply described by the numbers 
{riij + n'^j}- Here is a simple example: 







1 


1 


2 


2 


2 


3 




4 











1 




1 


2 




2 













1 


1 




1 


1 


2 


2 


2 




2 


2 


3 






A 













(4.5) 



This ring of tableaux, then, is isomorphic to the model R constructed above and we 
do not distinguish between them. Thus we specify a set of elementary couplings (i.e. a set 
of generators of R) as a set of elementary LR Tableaux. 



4.1- Example: the su{2) case 

The complete set of inequalities for su{2) variables {Ai,nii, 77-12} is simply 

Ai > ni2 nil > ni2 > (4.6) 

The other weights are fixed by 

fJ'i = nil + ni2 1^1 = Xi+ nil - "-12 (4.7) 

By inspection, the elementary solutions of this set of inequalities are 

(Ai,nii,ni2) = (1,0,1), (1,0,0), (0,1,0) (4.8) 

which correspond respectively to Ei,E2,E3 in (2.23). These correspond to the following 
LR tableaux: 



El 



Eo 



□ 



(4.9) 



It is also manifest that there are no linear relations between these couplings. The generating 
function is thus simply: 

1 



Qsu{2) _ 



{1-Ei){l-E2){l-Es) 



(4.10) 



(Recall once again our convention concerning grading and model variables: the LR tableaux 
in this generating function must be replaced by the corresponding monomials in the grading 
variables). 
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4-2. Example: multiple tensor products in the su{2) case 



Before we turn to more complicated algebras, it is interesting to reconsider the problem 
of multiple tensor products treated previously from the character method. Let us look for 
the multiplicity of the representation ( in the triple product A (8) (8) D C- In the first 
step, the LR rule applies as before: with rin + 77-12 = //i, we have Ai > 77-12. After the first 
product, we re-apply the LR rule with now Ai replaced by Ai +nii — ni2 and Uij replaced 
by rriij with 777-11 + "712 = i^i- The LR gives Ai + nn — ni2 > mi2- The two inequalities 
that defines the su{2) LR basis for the quadruple product are then 

Ai > ni2 Ai + nil — '"-12 > "712 ti^- > niij > (4-11) 

The elementary solutions are then, in the order: name of the coupling, corresponding 
Dynkin labels and the 5-vector (Ai, Tin, 7112, ttih, 77112),: 

El : (1)®(1)® (0) D (0) (1,0,1,0,0) 

£^2 : (1) (0) (1) D (0) (1,0,0,0,1) 

Eg : (1) ® (0) (0) D (1) (1,0,0,0,0) 

(4.12) 

E4 : (0) ® (1) ® (1) D (0) (0,1,0,0,1) 

E5 : (0) ® (1) ® (0) D (1) (0,1,0,0,0) 

Ee : (0) (g) (0) (g) (1) D (1) (0,0,0,1,0) 

The linear relation, whose existence was signalled by the character method, is 

E3E4 = E2E5 : (1, 1, 0, 0, 1), ^ EiEe : (1, 0, 1, 1, 0) (4.13) 

Choosing to forbid the product E3E4, the generating function can be written in the form 

1 — EsE4 



G = 



(1 - Ei){l - E2){1 - E^){1 - E4)(l - E^){1 - Eq) 



n t4 



+ 



l-Ei \ \l-E^ 1-E. 

The latter form makes manifest the absence of E^E^. 

We could represent the elementary couplings in terms of tableaux, where the boxes 
with 1 refers to the ji tableau and those with 2 originates from the v tableau. (Warning: the 
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resulting tableaux describing the four-products are not necessarily LR tableaux.) Hence, 
nij gives the number of 1 in row j of the composed tableau while mik gives the number 
of 2 in row k. The elementary tableaux are 



El : 
E4: 



E' 



2 '■ 



E3: □ 



(4.15) 



From this representation, the relation reads 

E^E/i = E2E^ : 



7^ EiE^ 



6 : 



(4.16) 



It is worth pointing out here the relation between the multiplicities of the su{2) triple 
tensor products and sw(3) weight multiplicities in sw(3) irreducible representations. To 
make this relation explicit, we first consider void boxes to be filled with and then reshuffle 
all filling numbers by 1, i.e., 

J - P (4.17) 

The multiplicity of the product \ ® n ® f ^ C, is then seen to be the number of tableaux 
filled by Ai I's , /ii 2's, z/i 3's (the tableau has thus a total number of Ai + /Ui + z/i boxes) 
with the restriction that there should be C,i columns of 1 box. Since by construction, the 
numbers are strictly increasing in each column from top to bottom and non-decreasing in 
each row from left to right, these tableaux are nothing but sw(3) semistandard tableaux. 
The resulting tableau has sw(3) weight^ 



\lUJl + Hi{-U}i + UJ2) - 1^1(^2 



(4.18) 



and row lengths: (Ai +^1 + 1^1+ Ci)/2 and (Ai + fXi + i^i — Ci)/2. In terms of Dynkin 
labels, the number of such semistandard tableaux is exactly the multiplicity of the su{3) 
weight (Ai — /Ui, /Ui — z^i) in the representation of highest weight (Ci, (Ai + fXi + ui—(i)/2). 
For instance, the su{2) product (1) (1) (1) D (1) has multiplicity 2, corresponding to 
the products (forbidding E2E5) and EiE^, or equivalently, to the two tableaux 



(4.19) 



10 



In a su{N) semistandard tableau, a box filled by i has weight —u)i-\+0Ji, with ujq 



0). 
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and this is exactly the multiphcity of the su{3) weight (0,0) in the highest-weight repre- 
sentation (1,1), whose two semistandard tableaux are 



(4.20) 

Equivalently, since the multiplicity of the su{2) product A®/UCg>z^ D Cis the same as 
X (E) (E) (E) C ^ (0), by filling the boxes of C with 4's, one can relate the resulting 
multiplicity to a su(4) weight multiplicity in an irreducible representation of two equal 
rows. More precisely, the resulting tableaux have Ai I's, /ii 2's, ui 3's and Ci 4's, which 
corresponds to a weight 



XlUi + 111 {-uoi -I- u;2) + z^i (-0^2 + ws) - Cll^3 
and the corresponding irreducible representation has weight (0, n, 0, 0) with 

2n = Ai + //I + + ^1 
More generally, the multiplicity of the multiple A'" > 3 su(2) product 



(4.21) 



(4.22) 



N 



X^^ D (0) 



(4.23) 



is equal to the multiplicity of the su{N) weight A' in the highest-weight representation A, 
with 



TV TV 

A' = E >^ (-^-1 + A = (0, i J] Af ) , 0, • • . , 0) 



(4.24) 



(with ujQ = CON = 0). 



4-3. Example: the su{3) case 

Let us return to standard triple tensor products and turn to su{3). The LR conditions 
for su{3) are 

Ai > ni2 



A2 > ni3 
A2 + ni2 > ni3 + n23 

nil > n22 
nil + ni2 > n22 + ^23 



(4.25) 
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The other weights are given by 



Ail = nil + ni2 + ni3 - n22 - "23 

/i2 = ri22 + n23 

ui = Ai +nii - ni2 - n22 

z^2 = A2 + ni2 + n22 - nis - ^23 



(4.26) 



The elementary solutions of the set of inequalities (4.25) are again easily found by in- 
spection and they are given by the following set of numbers: (Ai,A2,nii,ni2,ni3, 77-22, ^23] 



El : 


(1,0) 


® (0,1) 


D 


(0,0) : 


(1,0,0,1,0,0,1) 


E2 '■ 


(1,0) 


(0, 0) 




(1,0) : 


(1,0,0,0,0,0,0) 


E3 : 


(0,0) 


(1,0) 




(1,0) : 


(0,0,1,0,0,0,0) 


E4: 


(0,1) 


(1,0) 




(0,0) : 


(0,1,0,0,1,0,0) 


E5: 


(0,1) 


(0,0) 




(0,1) : 


(0, 1,0,0,0,0,0) 


Eq : 


(0,0) 


(0,1) 




(0,1): 


(0,0,1,0,0,1,0) 


i?7 : 


(1,0) 


(1,0) 




(0,1): 


(1,0,0, 1,0,0,0) 


Es : 


(0,1) 


(0,1) 




(1,0): 


(0,1,1,0,0,0,1) 



(4.27) 



The corresponding tableaux are 



El : 



E. : 



E2: □, ^3: E, E, 



Ef\ : 



E, 



Es : 



(4.28) 



If we use the 7-component vector as a description of an elementary coupling, we see that 
there is one relation 

E1E3E5 = E7E8 : (1, 1, 1, 1, 0, 0, 1) (4.29) 



In section 7, we present a way to generate the elementary solutions of a system of Diophantine 
inequalities starting from the construction of a vector basis. This method appears to be avoid some 
of the computational complexity of the MacMahon method and in many cases the calculations 
can be done by hand. 
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This is confirmed by the construction of the corresponding LR tableaux: 

EiE^E^ : I ■ E ■ 

and 









1 






1 














1 






1 




2 







(4.30) 



(4.31) 



The generating function can thus be constructed by either forbidding EiE^E^ or EjEg. 
That yields directly the first or the third expression of (2.46) respectively.lii 



4-4- Example: the su{4) case 

To the su{4) LR known conditions are 





Ai 


> 


ni2 




A2 


> 


ni3 


A2 


+ ni2 


> 


ni3 




A3 


> 


ni4 


A3 


+ nis 


> 


ni4 


+ ni3 


+ n23 


> 


ni4 




nil 


> 


n22 


nil 


+ ni2 


> 


n22 


+ ni2 


+ ni3 


> 


n22 




n22 


> 


n33 


n22 


+ ^23 


> 


n33 



(4.32) 



The tensor- pro duct elementary couplings, that is, the elementary solutions to these 
inequalities are best written directly in terms of LR tableaux: 



Ai: 



A3:U: Bi: 



Ci:B, C2: 



C3 : 











, D'l-. 











B. : 



D'2: 



S, : 



D'3-- 



(4.33) 



12 



In this basis, E2E4EQ is an independent product. 



30 



together with 



El 







1 






1 


D2 : 




2 


, Ds: 








3 






2 





Ez '■ 





1 




2 


1 




3 





(4.34) 



The Dynkin-label transcription of the elementary couplings reads 



Ai : 


(0,0,0) 


® (0,0,1) 


D 


(0,0,1) 


D[: 


(0,1,0) 


® (1,0,0) D 


(0,0,1) 


A2: 


(0,0,1) 


(1,0,0) 




(0,0,0) 


D'2: 


(1,0,0) 


® (1,0,0) D 


(0,1,0) 


A3: 


(1,0,0) 


® (0,0,0) 


D 


(1,0,0) 




(1,0,0) 


® (0,1,0) D 


(0,0,1) 


Bi : 


(0,0,0) 


«) (0,1,0) 


D 


(0,1,0) 


Di : 


(0,1,0) 


(8) (0,0,1) D 


(1,0,0) 


B2 : 


(0,1,0) 


(8) (0,1,0) 


D 


(0,0,0) 


D2 : 


(0,0,1) 


(8) (0,0,1) D 


(0,1,0) 


B3: 


(0,1,0) 


® (0,0,0) 


D 


(0,1,0) 


D3: 


(0,0,1) 


® (0,1,0) D 


(1,0,0) 


Ci : 


(0,0,0) 


(1,0,0) 


D 


(1,0,0) 


El : 


(1,0,1) 


(0,1,0) D 


(0,1,0) 


C2: 


(1,0,0) 


® (0,0, 1) 


D 


(0,0,0) 


E2 : 


(0,1,0) 


(0,1,0) D 


(1,0,1) 


C3: 


(0,0,1) 


® (0,0,0) 


D 


(0,0,1) 


E3 : 


(0,1,0) 


(8) (1,0,1) D 


(0,1,0) 



(4.35) 



For stt(4), there is a large number of linear relations: in fact there are 15 relations [12,5] 

EiEj = BkDkDf, 



CiEi 
DiEi = CjBf.Dk 



DjDf, — BiCjCk 
DiEi = BjDjCk 



(4.36) 



with i, j, k a cyclic permutation of 1,2, 3. 



To construct the generating function, we need to select forbidden couplings. It turns 
out that when there are more that one relation, complications may arise. We must ensure 
that the selected forbidden couplings are complete, which means that no further (usually 
higher-order) relations are required for a unique decomposition of a given coupling. How do 
we select a set of complete compatible forbidden couplings? A technique that is tailor-made 
for dealing with problems of that type is that of Grobner bases. This will be introduced in 
the next section. At this point, we simply indicate a complete choice of forbidden couplings, 
namely {EiEj, D'^Ei, DiEi, DjO^, DjDi}. This yields then a model for the generating 
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function, which then reads [12,5] : 

3 

+ D2D2D3E1 + D1D1D2D3 + E3E3D1D2 + D[D[DiE3 (4.37) 

+ D'^Es&^EsD'^ + E2E2D[D'3 + E2DiE2DiD[ + £;2A3^2^3^3 
+ D1D3E2D1D3E2 + D2D'^D2D2Ei + D2D'^E3D2D'2E3). 

where 

Mi = {l-Mi)-\ (4.38) 

4.. 5. A remark on the reduction of grading variables and higher multiplicities 

In this section, we would hke to stress the fact that a generating function built from a 
complete basis (or a complete description of tensor products) has 'multiplicity coefficients' 
all either or 1: the 'coupling' either exists or not. Higher multiplicities can be gener- 
ated only after the number of grading variables has been reduced (that is, some grading 
variables have been set equal to 1). To be explicit, suppose that we consider the gener- 
ating function for the su{N) tensor products starting from a LR description. This is a 
system of inequalities for the variables {Aj, rijk} with z, j = 1, • • • , A?" — 1 and j <k =< N. 
Introducing the grading variables {Li, Njj.}, the generating function reads: 

G = J2C{x.,n,.}L}^N;^'' (4.39) 

and the coefficient C^x^^njk} is either 1, if the solution exits, or 0, if the solution does 
not exist. In other words, when expressed in terms of a complete set of variables, the 
multiplicity is or 1. Larger values for the multiplicity can only result from the reduction 
in the number of grading variables. Typically, we consider the reduction from the set 
{Xi,njk} to the set {A^, /i^, z/^}. The multiplicity, which now has the interpretation of a 
tensor-product coefficient, is then no longer trivially or 1 (except for su{2) where this 
changes of variables does not induce a reduction). 
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5. Diophantine inequalities: elementary couplings, relations and Grobner bases 



We will introduce the idea of the Grobner basis via a simple example. Suppose 
we know that R = Q[x, y, z,t\/ 1 where / = {xy — t, zy — t), with an grading given by 
(1, 0), (0, 1), (1, 0) and (1, 1) for x, y, z and t, is a model for a generating function. Writing 
X = x-{-I and similarly for the other variables, we have in R that xy = t and zy = t. These 
two expressions give two re- write rules : xy ^ t and zy ^ t. These rules can be used to 
simplify any monomial. The aim is to find a re-write rule which, when iterated, produces 
unique representatives for the classes of /. If this is the case, then a vector space basis of 
R would consist of terms of the form m + I with m a monomial which is not divisible by 
any of the left-hand sides of the rewrite rules. 

In the example above, if we had 'good' rewrite rules then a basis for R would be 
represented by monomials not containing xy or zy, i.e. monomials of the form either y°'t^ 
or x°'z^t^. The generating function which counts there monomials, taking into account the 
grading and potential over-counting, is: 



The exponent of A carries the first grading index and B the second.EJ 

However this generating function is not correct. It contains the term 2A^B corre- 
sponding to the 2 monomials xt and zt. But the polynomial z{xy — t) — x{zy — t) = xt — zt 
is also in / and hence in R we have xi = zi and so the space of grade (2, 1) has dimension 
1 rather than 2. This problem can also be seen as a problem with the re- write rules. If we 
start with xyz then we can use the first re-write rule: xyz ^— tz or the second: xyz ^— xt. 
We cannot apply any further re- write rules and so this set of re-write rules does not pro- 
duce a unique representative. The solution is to include the rule xt i— > zt. This gives a set 
of 3 rules: xy ^ t, zy ^ t and xt ^ zt. It turns out that this is a 'good' set and so a basis 

For an elementary introduction to Grobner bases, see for instance [13]. 

In more details: with the specified grading, y°'t^ corresponds to the term B°'{AB)^ which is 
generated by [(1 — AB){\ — B)]^^ while x°'z^t'^ corresponds to A"" A^{ABY which is generated by 
[(1 — AB){1 — AY]^^ . Since the constant term would be counted twice if we simply add these two 
pieces, we multiply the first one by B. 




(5.1) 
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for R is given by (the classes of) monomials of the form y°'t^, x^'z^ and z^'t^ which gives 
the generating function: 

' + 7tA^ + (5-2) 



{1- AB){1- B) (1-A)2 {1-A){1-AB) 

The set of 'good' generators, xy — t,zy — t, xt — zt we have found for / is known as a 
Grobner basis [14]. 

The general procedure for constructing a Grobner basis given a set of generating 
polynomials is as follows. First choose a term ordering, which is an ordering on monomials 
with the property that any chain mi > m2 > . . . has finite length. For example we can 
order the variables hy x>y>z>t and then order all monomials by the corresponding 
lexicographic (dictionary) order, for example: x'^y > xyz > y^. For each generator of our 
ideal /, select the monomial which is highest with respect to the given term ordering. This 
is then the term which appears on the left of the re-write rule. The lexicographic ordering 
gives the first two re-write rules of our example: xy ^ t and zy i— > t. Next, for each pair 
of leading terms find the lowest common multiple and simplify it in the two possible ways. 
In this case there is only one pair of leading terms and the lowest common multiple is xyz 
which simplifies to xt and yt. Continue to apply the re-write rules until the terms do not 
simplify further. If the resulting pair of terms are the same, then proceed to the next pair 
of leading terms, otherwise add a new re- write rule. In this case we add xt yt. Proceed 
until no pair of leading terms gives a new rule. This is the case for the rules we now have. 
For example the two rules xy ^ t and xt ^ zt appears to give a new rule by simplifying 
xyt to both t^ and yzt. However the second term can be further reduced to and so 
no new rule is required. Improvements on this basic algorithm, known as the Buchberger 
[14] algorithm, mean that it is now feasible to find Grobner bases for quite large sets of 
generating polynomials.^ 

Although it is not clear from this example, the technique of Grobner bases is a very 
versatile tool for performing explicit calculations. We end this section with an illustrative 
example relevant to our discussion of tensor-product generating functions. 

Consider a set of linear Diophantine equations: 

Ma = 0, a > (5.3) 



The web pages of the computer-algebra information network at |http:/ /cand. can.nl/CAIN 
contain information about many of the programs currently available. 
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with M an integer matrix and a a vector of non- negative integers. We would like to 
construct a generating function for the solutions to this set of equations: 



E 



X 



(5.4) 




A non-trivial example is given by the Diophantine equations that describe a 3 x 3 magic 
square: 

fa b c\ 

(5.5) 

with non-negative entries and equal row and column sums. The magic square condition 
(the sum of each row and each column is the same, say equal to t) gives the following set 
of equations: 

a + b + c = t 
d + e + f = t 
g + h + i — t 
a + d + g = t 
b+e+h=t 
c + f + i = t 

With a standing for the column vector with entries (a, b, c, d, e, /, g, h, i, t), the matrix M 

reads 

y^l 1 1 -IN^ 
1110 0-1 
111-1 
10 10 10 0-1 
10 10 10-1 
Vo 1 1 1 -1/ 



(5.6) 



M 



(5.7) 



It can be show than any solution of Ma = is a linear combination (with non-negative 
coefficients!) of a finite number of basic solutions cki, • • • , cug (see for example [4]). Moreover 
there is a straightforward algorithm for finding this basic set of solutions [15]. In this case 
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we have the following set: 

ai = (0,0,1,0,1,0,1,0,0,1) 
^2= (0,1,0,0,0,1,1,0,0,1) 

a3 = (0,0,1,1,0,0,0,1,0,1) 

(5.8) 

a4 = (1,0,0,0,0,1,0,1,0,1) 
as = (0,1,0,1,0,0,0,0,1,1) 
ae= (1,0,0,0,1,0,0,0,1,1) 

We shall use A, S, . . . , T to denote the "grading variables" of this example so that the 
exponent of A carries the value of a and so on. A model for the generating function is given 
by the subring S of Q[^, B, C, D, E, F, G, H, /, T] generated by monomials corresponding 
to the 6 elementary solutions. 

El = CEGT, E2 = BFGT, E3 = CDHT, 

E^ = AFHT, E5 = EDIT, Eq = AEIT 

The monomials in S correspond to magic squares. For example EIE4EQ = A^C'^E^FG'^HIT^ G 
S corresponds to a square with row and column sums equal to 4: 

2 2\ 

3 1. (5.10) 
2 11/ 

Note that in this example it is convenient to construct our model as a subring of the 
ring of grading variables. Thus each "elementary coupling" Ei is actually equal to the 
corresponding monomial in the grading variables. 

However, there are relations between these generators and so it is not immediately 
clear how to construct the Poincare series for S. What we require is an isomorphism of S 
with R = Q[ei, . . . , ee]// such that 1— > z = 1, . . . , 6 and such that we have a Grobner 
basis of the ideal / (the 'ideal of relations'). 

Fortunately, such an isomorphism is easily constructed using Grobner-basis meth- 
ods. Introduce the ring Q[A, S, C, D, i?, F, G, i/, /, T, . . . , Eq] with the lexicographic 
ordering 

A>B>C>D>E>F>G>H>I>T>ei>...>eQ (5.11) 

Let J be the ideal generated by Ei — ei, . . . , Eq — cq. This is not necessarily a Grobner 
basis with respect to this term ordering. Let G be the Grobner basis for J with the given 
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ordering. Then it can be shown [14] that G fl Q[ei, . . . , ee] is a Grobner basis for the ideal 
of relations / which we require. In this case G is quite large, but its intersection with 
Q[ei, . . . , ee] is 616465 — 626365. The corresponding relation in R is E1E4E5 — E2E^Eq and 
it is straightforward to verify that these two terms do indeed give the same magic square, 
so that indeed we have found a relation between the generators of R. The Poincare series 
for Q[ei, . . . , ce]// is easily computed: 



1 f I 



{1 - E2){1 - Es){l - Ee) \{1 - Ei){l - E4) 



_l E5 ^ E4E5 



{1-E^){1-E5) {1-E4){1-E5) 



(5.12) 



If we do not need to keep track of the row and column sums we can set T = 1 in this 
generating function. Similarly, if we simply want a generating function for the number of 
magic squares with given row and column sums, we can set all variables except T equal to 
1, which yields: 

1 _|_ y _|_ ^2 

Z ■ (5.13) 

(1-T)5 ^ ^ 

It is perhaps worth pointing out that the model we constructed above is still a model for 
(5.13) if we change the grading by keeping just the last component of the grading vector. 
Now, however, the grading variable is simply T and so we cannot identify the grading and 
model variables. Thus to obtain (5.13) from (5.12) we once again must use our convention 
that model variables are replaced by the corresponding grading variables. In this case this 
corresponds precisely to setting all the variables except T equal to 1. 

We have made use of the techniques described in this section to verify the results we 
have quoted previously. 



6. The Berenstein-Zelevinsky basis 

6.1. Generalities 

The presentation of the stt(2, 3, 4) cases has made clear the usefulness of a basis, or 
more precisely, a re-expression of the tensor-product calculation in terms of Diophantine 
equations. The Littlewood-Richardson algorithm yields a set of such equations only for 
su{N). Fortunately, Berenstein and Zelevinsky [16] have expressed the solution of the 
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multiplicity of a given tensor product as a counting problem for the number of integral 
points in a convex polytope. For a given algebra, the polytope is formulated in terms of a 
characteristic set of inequalities. These inequalities can be interpreted as defining a tensor- 
product basis. For su{N), this basis reduces to the LR set of inequalities. For the other 
classical algebras but sp{4), the proposed basis is a conjecture. Hence, in order to analyse 
our next example, sp{4), we first need to present the corresponding Berenstein-Zelevinsky 
inequalities. 



6.2. The BZ sp{4) basis 

The combinatorial description of tensor products for sp{4) is not as simple as in the 
standard LR product must be supplemented by a division operation and 
modification rules [17]. Various natural trials for the search of a linear basis failed (and 
for instance, a stretched product making manifest the relations between the elementary 
couplings cannot be defined). 

Given the BZ set of inequalities, the natural way to proceed, as just mentioned, is to 
interpret these as the appropriate basis for the description of the tensor products. These 
inequalities are as follows 

Ai > Pl2 

A2 > ri/2 

A2 > ri/2 + qi2 -pi2 
A2 > r2/2 + qi2 -Pi2 
/Ui > qi2 

/Ui > qi2 +ri-r2 
fJ-i > P12 + ri - r2 
IJ'2 > r2/2 

vi = r2 - ri - 2pi2 + Ai + ;Ui 

T^2 = P12 - qi2 - r2 + A2 + /X2 



(6.1) 



Our notation is different from that used in [16]. The relation is ri = mi, r2 = m2, P12 
mi2, qi2 = rnX^. 

38 



The sp{4) tensor product coefficient Afx^i, is thus given by the number of solutions of the 
above system with ri, r2 G 2 Z_|_ et pi2, qi2 G Z_|_. 

A proper set of variables for a complete description of a particular tensor-product 
coupling is thus {Ai, A2, /^i, /i2, ^'i, ?'2,Pi2, q'i2}- We give the list of elementary couplings, 
adding to each coupling the corresponding four- vector [ri, r2, P12, ^12]: 



Ai : 


(0,0) (g) 


(1,0) 


D 


(1,0) 


[0,0,0,0] 


A2 : 


(1,0) o 


(0, 0) 


D 


(1,0) 


[0,0,0,0] 


^3 : 


(1,0) ® 


(1,0) 




(0, 0) 


[0,0,1,1] 


Bi : 


(0,0) ® 


(0,1) 


D 


(0,1) 


[0,0,0,0] 


B2: 


(0,1) ® 


(0,0) 


D 


(0,1) 


[0,0,0,0] 


B3: 


(0,1) ® 


(0,1) 


D 


(0,0) 


[2,2,0,0] 


Ci : 


(0,1) ® 


(1,0) 




(1,0) 


[0,0,0,1] 


C2: 


(1,0) ® 


(0,1) 


D 


(1,0) 


[0,2,1,0] 


C3: 


(1,0) ® 


(1,0) 




(0,1) 


[0,0,1,0] 


Di 


(2,0) ® 


(0,1) 




(0,1) 


[0,2,2,0] 


D2 


(0,1) ® 


(2,0) 


D 


(0,1) 


[2,0,0,0] 


Ds 


(0,1) ® 


(0,1) 




(2,0) 


[0,2,0,0] 



(6.2) 



The unspecified linear relations mentioned in (2.59) can now be obtained. To find 
those products that are really equal in the present basis, we only need to compare their 
corresponding sets of four-vectors [ri, r2,pi2, ^12] (which are additive in products of cou- 
plings). We thus find for instance that 



C1C2 = AsDs : [0, 2, 1, 1] ^ Ai A2S3 : [2, 2, 0, 0] 



(6.3) 



Proceeding in this way for the other cases, we find the following complete list of relations: 
C1C2 = AsDs, C2C3 = AiDi C3C1 = A1A3B2 

DiD2 = BsCi D2D3 = AIB2B3 DiDs = B2CI (6.4) 

CiDi = ASB2C2 C2D2 = A1S3C3 C^Ds = A1B2C2 

The use of the BZ basis to find the elementary couplings and their relations seems to 
be novel. 



"""^ Our construction provides an indirect proof of the validity of the BZ basis in that from it we 
recover the result of [10] derived from the character method. 
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A possible choice of forbidden products is the one given in [10]: 



(6.5) 



with z, J = 1, 2, 3 and i ^ j- It leads to the generating function: 



+C2C2D3 + D^CsD^ + CsCsD^ 



) 



(6.6) 



Of course, by modifying the ordering in the Grobner basis, we can get other choices of 
forbidden couplings. Here is another set of forbidden couplings that can be obtained: 
{DiDj,CiDi, AiDi, A'iD^, ^1^352}- The corresponding generating function reads 



QSp{i) ^ B1B2BS 



(6.7) 



+D1D1A2A3C2C3 + D2D2AiA2A3CiC3{l - A1A3B2) 



These two generating functions are equivalent when rewritten in terms of the grading 
variables, that is, in terms of Dynkin labels. However, they originate from two distinct 
models. The second one turns out to be well adapted to the fusion extension. 



7. A vector-basis approach to the construction of generating functions 

In this section, we present a simple and systematic way of generating by hand all the 
elementary solutions of a set of linear homogeneous inequalities starting from the well- 
known construction of a vector basis. The first step amounts to reformulate the system of 
inequalities in terms of equalities. We then look for the elementary independent solutions 
by relaxing the positivity requirement. In other words, we construct the vector basis. In a 
final step, we find the minimal linear combinations of these vector-basis elements that yield 
positive solutions. This will also provide a further illustration of the MacMahon's projec- 
tion technique. But actually, the outcome of this projection is the desired tensor-product 
generating function. Hence, this approach turns out to be a new way of constructing the 

40 



generating functions. Ill This method wiU play a key role in our subsequent description of 
fusion-rule generating functions. 

In the first subsection, we present the equality version of the LR system of inequalities 
which holds for su{N) . This happens to lead to a remarkable graphical description of su{N) 
tensor products that has been found by Berenstein and Zelevinsky [18]. In the following 
section, we illustrate the procedure by working out the vector-basis argument for 5^(3). 
A more general discussion of this technique is presented in the next section. We finish 
with an analysis of the sp{4) case, for which we provide a novel graphical representation 
of tensor products. 



7.1. su{N) tensor-product basis: from inequalities to equalities: graphical representations 
as BZ triangles for su{N) 

Consider the direct transformation of the LR inequalities to equalities by introducing 
an appropriate number of new non-negative integer variables. Consider first the su{2) 
case, for which there is a single inequality: Ai > 7112. We transform this into an equality 
by introducing the positive integer a defined by 

Ai = ni2 + a (7.1) 

The expression for ui becomes then ui = Xi + nn — 77,12 = a + uu. Since /ii = nn + 71-12, 
we are led naturally to a triangle representation of the tensor product: 

X^LiDu ^ ^ (7.2) 

ni2 nil 

We read off the Dynkin label of the A representation from the sum of the two integers 
that form the left side of the triangle, that of the n representation from the bottom of the 
triangle and the ui label is the sum of the two integers that form the right side. A more 
uniform notation amounts to setting a = mi2 and Tin = /12, in terms of which the triangle 
looks quite symmetrical: 

"-12 n2 



This method, referred to as being novel, is probably well-known but we could not trace it 
precisely in the literature. There are implicit remarks in that direction in the first reference of [4]. 
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with 

Ai = mi2 + ni2 fJ,i = ni2 + ^12 1^1 = mi2 + I12 (7.4) 

These numbers mi2 and I12 plays the role of ni2 in the permuted versions of the tensor 
product. The triangle combinatorial reformulation of the tensor product problem is as 
follows: the number of triangles that can be formed from nonnegative integers 77-12, "7-12 
and li2 that add up to the Dynkin labels of the representations under study according to 
the above relations gives the multiplicity of the triple coupling A (g) D i', or equivalently, 
the multiplicity of the scalar representation in the product \ D (0) (since for stt(2), 

V* = u). 

For su{3) the situation is somewhat more complicated. The transformation of the 
inequalities (4.25) into equalities takes the form 

Ai = ni2 + a 

A2 = ni3 + 6 

A2 + ni2 = 77-13 +n23+c (7.5) 

nil = n22 + d 
nil + ni2 = n22 + ?^23 + e 

The expression for the other weights becomes 

A*! = ni3 + e 112 = n22 + n23 

(7.6) 

vi = a + d V2 = n22 + c 

Since there are two expressions for both nn and A2, there follows the compatibility rela- 
tions: 

ni2 + d = 7123 + e 7123 + c = 6 + ni2 {1.1) 
By adding these two relations, we find: 

c + d=6 + e (7.8) 

Again we are led naturally to a triangle representation: with C = this reads 



a 

ni2 d 
b c 
ni3 e 7123 n22 



(7.9) 
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We read the Dynkin labels from the sides of the triangles, from Ai to C2 in an anti-clockwise 
rotation starting from the top of the triangle, exactly as for su(2), except that here there 
are two labels on each sides. Notice that the compatibility conditions amounts to the 
equality of the sums of the extremal points of the three pairs of opposite sides of the 
hexagon obtained by dropping the three corners of the triangle. 

Again a more symmetrical notation calls for itself: 

a = mi3 b = m23 c = mi2 d = I23 e = I12 n22 = I13 (7-10) 



in terms of which the triangle reads 



mi3 
ni2 hs 
"7,23 rni2 
ni3 l\2 n23 h3 



with labels fixed by: 



(7.11) 



Ai = mi3 + ni2 A2 = m23 + ni3 

fJ'i = ni3 + li2 /i2 = n23 + h3 (7-12) 
Ci = li3 + mi2 C2 = I23 + mi3 
The hexagon conditions read: 

ni2 + m23 = ^23 + mi2, 

^12 + "7,23 = ^23 + "7,12, (7-13) 
ll2 + "23 = ^23 + "12- 

In terms of triangles, the problem of finding the multiplicity of the su{3) tensor product 
A (8) (g) D boils down to enumerating the number of triangles made with nonnegative 
integers that form a bipartition of the Dynkin labels and that satisfy the above three 
hexagon relations. 

For stt(4) the BZ triangle is obtained in a similar way. One first transforms the 
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inequalities (4.32) into equalities by introducing positive integers a,b, - ■ • , k: 



Ai 


= ni2 + a 


A2 


= ni3 + 6 


A2 + ni2 


= ni3 + n23 + c 


A3 


= ni4 + d 


A3 + ni3 


= ni4 + 7124 + e 


A3 + ni3 + n23 


= ni4 + n24 + ^34 + / 


fill 


— ''■22 y 


nil + ni2 


= n22 + "-23 + h 


nil +ni2 + nis 


= n22 + n23 + n24 + i 


n22 


= ?^33 + j 


n22 + n23 


= ^33 + 77,34 + K 



(7.14) 



The Dynkin labels of fj, and C = ^* becomes 

A*! = ni4 + i, fX2= n2A + M3 = "-33 + "-34 

(7.15) 

Ci =^33 + /, C2 = c + J, C^ = a + g 
Since every Dynkin label is the sum of two positive integers, a triangle representation is 
again natural: 

a 

ni2 g 

^ (7.16) 
y V 

d f 

ni4 i s t n34 n33 

However the position of the integers specifying the labels 2 is ambiguous at this point: x 
can be either b or nis {y being the other one), similarly {u,v) is related to the doublet 
(j, c) and (s,t) to (7124, «;)• Moreover not all the needed integers appear in this triangle: 
e, h, 77-23 are missing. In order to take into account the various compatibility relations, it is 
natural to insert the three remaining points in the center of the big triangle forming then 
three hexagons: 



(7.17) 

• • • • ^ ^ 

• • • 

• ••••• 
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Quite remarkably, the different conditions are simply the equality of the opposite sides of 
every hexagon, exactly as for su{3), and these completely fix the position of every integers 
in the triangle. For instance, replacing A2 = ni3 + b into A2 + ni2 = ni3 + 77-23 + c yields 

ni2 + b = n23 + c (7.18) 

For this to be an hexagon relation, b and c must belong to the same hexagon and since 
they are on opposite sides of the triangles, they must belong to the first hexagon: hence 
X — b and u = c. That also fixes the position of ?i23 inside the big triangle. Proceeding in 
this way with the other constraints, we end up with the following representation: 

a 

ni2 g 

^ , . (7.19) 

ni3 h n23 J 

d e f 

ni4 i n24 « ri34 7733 

Finally, a more symmetrical form is obtained by redefining the name of the positive 
integers a, - ■ ■ ,k and 77-33 as follows: 

"^14 
ni2 ^34 
77724 "713 
^13 I23 ^23 I2A 
77134 77123 "^12 

7li4 li2 7124 ^13 "-34 ^14 



(7.20) 



They are related to the Dynkin labels by 

Ai = mi4 + ni2 A2 = 77724 + "■13 A3 = 7/734 + 7li4 

Atl = niA + I12 A*2 = "24 + ^13 A*3 = "-34 + ^14 (7-21) 

Ci = liA + rni2 C2 = ^24 + rni3 C,3 = I3A + 77214 
and the hexagon relations read: 

"-12 + TO24 = rni3 + 7123 "-12 + ^34 = ^23 + "-23 "^24 + ^23 = ^34 + TO13 

?^13 + ^23 = h2 + "-24 "-13 + "^34 = "-24 + "^23 "^34 + ^12 = ^23 + "^23 (7.22) 

^24 + ?^23 = ^13 + ^^34 "-23 + "^23 = "^12 + "-34 ^13 + "^23 = ^24 + "^12 

The su{N) generalisation is obvious; the triangle is built out of {N — 1){N — 2)/2 
hexagons and three corner points. 
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Here is the rationale for the labeUing nij,mij, lij from the triangle point of view [18]. 
If Ci are orthonormal vectors in R^, then the positive roots of su{N) can be represented 
in the form Cj — e-,-, I < i < j < N. The triangle encodes three sums of positive roots: 

x + fi-C 

The hexagon relations are simply the consistency conditions for these three expansions. 
Clearly, the variables riij that appear in the above relations are exactly the that appear 
in the LR tableaux for the product X® fi D = v. 

7.2. Prom a vector basis to the generating function: the su{3) case 

Given the transcription of inequalities into equalities, we can easily extract the cor- 
responding basis vectors. This is the starting point of a new method for constructing the 
tensor-product generating functions. To keep things concrete, we focus on the su{3) case. 
The goal is to first get a vector basis and then to project it to get the elementary couplings. 
The generating function is a direct result of this procedure. 

The equality version of the LR inequalities have already been presented in the previous 
subsection: these are (7.12) and (7.13); they underlie the construction of the BZ triangle 
(7.11). The last hexagon condition of (7.13) is the difference of the previous two so it is 
not an independent relations. We thus have a total of 15 variables: Xi, ■■■ ,1/2,112, ■■■ , ^^23 
and 8 equations. The number of independent variables is thus 7. These will be chosen to 
be mis, 'm23,li3,l23,ni2,n'is, 7123- The dependent variables are fixed as follows: 



Ai 


= mi3 + 77,12 


A2 


= 17123 + ni3 


1^1 


= ni3 + ni2 + I23 - n23 




= n23 + ^13 


Ci 


= ni2 + m23 + ^13 - "-23 


C2 


= ^23 + mi3 


h2 


= ni2 + ^23 - "-23 


mi2 


= ni2 + 777,23 - ^^23 



(7.24) 



ej) , 

i<j 
i<j 

77y (cj — ej) , 



i<j 



(7.23) 
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We now look for the elementary solutions of this system (without invoking the constraint 
that all the above dependent variables should be necessarily positive). The sought basis 
vectors are obtained by setting one of the variable mi3, • • • , 71-23 to 1 and all other set equal 
to zero. This produces (in order) the triangles E2., E5, Eq, E^, E^, E4 and Zi displayed 
below: 

E2:(l,0)(0,0)(0,l) £;3:(0,0)(1,0)(0,1) 

1 
1 



1 

E^ : (0,1)(1,0)(0,0) E^ : (0, 1)(0, 0)(1, 0) Eq : (0, 0)(0, 1)(1, 0) 





11 

1000 0000 0001 

Er : (1,0)(1,0)(1,0) Zi : (0, 0)(-l, 1)(-1, 0) 


10 

1 0-1 
1 - 1 1 



(7.25) 



These are all genuine BZ triangles except for Zi which has some negative entries. However, 
at this level, there are no relations between these elementary solutions (the basis vectors 
are independent), hence the decomposition of any solution in terms of these 7 basic ones 
is unique. All solutions are then freely generated from the following function: 

^ " (1 - E2)il - Es){l - E4)(l - E5){1 - Ee){l - Er){l - Zi) ^^"^^^ 

To recover the generating function for all tensor products from the above expression, we 
need to project out terms that lead to triangles with negative entries. To achieve this, 
we introduce the grading variables associated to the above couplings (compare the above 
triangles with the general form given in (7.11)): 

E2 : Mi3, Es : L12L23 E4 : A^i3 

E5 : M12M23 Eq : Lis Er : L12M12N12 (7.27) 

Zi : L-'M-iiV23 

Our generating function follows from the projection of the above function G, re-expressed 
in terms of the grading variables, to positive powers of L12 and M12. Equivalently, one can 
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rescale L12 by x and M12 by y and project to positive powers of x and y and set a; = y = 1 
in the result. This is equivalent to the rescaling 

Es^xEs E^^yE^ Er ^ xyEr Z^^x'^y'^Zi (7.28) 

We are thus led to consider 

nnG{E2,xEsr--,x-^y-^Zi) (7.29) 
Keeping only those terms which depend explicitly upon x or we have then 



X y 1 



>> (1 - xEs){l - yE^){l - xyE7){l - x-^y-^Zi) 

1 



1 / 1 x-^y-^Zi 



+ 



{l-xE3){l-yE5){l-ErZi) \l - xyE^ 1 - V'^i 
No more work is need for the first term. For the second one, we have 

X y x~^y~^Zi f yE^ 1 



>> {1 - xE3){l - E7Zi){l - x-^ZiE^) \l-yE5 l-x-^y-^Zi 

X X~^EkZi 

O 



> (1 - E^){1 - E7Zi){l - xEs){l - x-^ZiE^) 
X x-^E^Zi ( xEz 1 



(7.30) 



(7.31) 



> {l-E^){\-E7Z^){\-EzE^Zr) \l-xE3 l-x-^ZiE^ 
_ E^E^Zi 
~ (1 - E^){1 - ErZi)il - EsE^ZiXl - E^) 

We then introduce the following two new elementary couplings 

El = E^Zi Eg = E^E^Zi (7.32) 
Collecting the two terms resulting from the projection, we end up with 

G*"^'^ = (jl (1 - ErEs) (7.33) 

which is indeed the su{3) tensor-product generating function. 

As a side remark, we indicate how the various triangles associated to a given triple 
product are related to each others. For this question, we consider the three weights to be 
fixed. The system of equations is now nonhomogeneous and solutions are given by a linear 
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combination of a particular solution and the sum of all homogeneous solutions, that is, 
solutions of the system with all Dynkin labels set equal to zero: 



=mi3 + ni2 
=m23 + ni3 
=ni3 + ni2 + hs - 1^23 

=77,23 + ^13 

=ni2 + m23 + ^13 - n23 

=^23 + mi3 

=ni2 + ^23 - ^23 - I12 

=ni2 + m23 - n23 - mi2 
These solutions are given by n A with n e Z and 

1 

A= -1 -1 
-1 - 1 

1 - 1 - 1 1 



(7.34) 



(7.35) 



Hence, given an allowed triangle, all other triangles related to the same tensor product can 
be obtained by adding or subtracting a number of times A while ensuring that all triangle 
entries are positive [19]. 

The great advantage of the vector-basis procedure is that the basis vectors are very 
easily obtained and their number grows slowly with N for su{N) (their number being 
{N — 1)(A'^ + 4)/2) as compared to the number of tensor-product elementary couplings: 



algebra 


# elem. coupl. 


dim. vect 


5-^(2) 


3 


3 


su{3) 


8 


7 


su{A) 


18 


12 


su{5) 


45 


18 


su{6) 


138 


25 


su{7) 


526 


33 



(7.36) 
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1.3. General aspects of the vector-basis construction 



In general, of course, the fundamental solutions to the linear system may have non- 
integral values of the variables. However the corresponding terms in the generating func- 
tion can be eliminated by rationalising all the denominator terms and then keeping only 
those terms in the numerator that have integral exponents. This suggests the following 
modification of MacMahon's algorithm. 

Consider the system of equations 

Ma; = 0, xeN'' (7.37) 

where M is a matrix of rank s. We thus have k variables and s relations between them. 
The dimension of the vector basis is thus k — s. We will denote the independent (free) 
variables as Xi, i = 1, ■ ■ ■ ,k — s and the remaining ones as Xj, j = 1, ■ ■ ■ , s. To find a 
generating function for the solutions of this system: 

1. First construct a basis in Q'^ for the solutions of Mx = by setting Xi = 1 with all 
other Xj zero (j = 1, ■ ■ ■ , — s, j 7^ i). Denote by x^^"* the value of the dependent 
variable Xj evaluated aX xi — 1 with all other Xi zero. The basis then reads 

ei = (l,0,0...,0;{£5')}), 
e2 = (0,l,0...,0;{xf }), 

(7.38) 

efe_, = (0,0,0,...,l;{£f-^)}) 

By construction, the e^'s are linearly independent. However notice that in general the 
x^j^ might be rational. 

2. From the form of the e^'s, it follows that any solution to (7.37) can be written as 

CiSi with Ci non- negative integers. In particular this means that every solution to 
(7.37) corresponds to a term in the generating function: 

^^^^ = ii-x^^){i-x^^)...{i-x^s) ^^-^^^ 

where Xi, . . . , Xj^ are grading variables. 

3. G{x) may contain negative or fractional exponents due to the occurrence of x^^^ in the 
exponents. These are eliminated by first using MacMahon's algorithm to eliminate 
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any negative exponents and then rationalizing denominators and keeping only terms 
with integral exponents in the numerators. 

The result is the generating function for the solutions to (7.37). This algorithm, 
however, does not seem to be optimal in all case. Consider the following example: 

3x-2y + z-3t = 

(7.40) 

2x + y-2z-t^0 

with x,y,z,t e N. Using linear algebra, we find two basic solutions: (|,|,1,0) and 
(|, — 1,0, 1). This means that the initial generating function is 



(1 - X3/7F8/7^)(^i _ X5/7F-3/7T) 



from which we must eliminate negative and fractional exponents. This is somewhat lengthy. 
However, the calculation can be shortened by some observations. First note that we can 
take as the fundamental solutions: (^^y^l^O) and (^5— y^O, |) since once again any 
solution to (7.40) is a non-negative linear combination of these two solutions. This give as 
the initial generating function: 

(7.41) 



(1 - X3/56y 1/7^1/8) (1 _ x5/21y-l/7yl/3) 

Keeping the positive exponents in Y yields: 

1 



(1 - X7/24^1/8rl/3)(l _ x3/56y 1/7^1/8) 

and keeping integral exponents in T and Y gives: 

1 

(1 - X7/8Z3/8T)(1 - X3/8rZ7/8) " 



(7.42) 



(7.43) 



Rationalizing these denominators and keeping integral exponents in the numerator yields 
the final generating function: 

(1 + X^Z^TY^ + X'^Z^T^Y^ + X^Z^T^Y+X^Z^T^Y'^+ 

(7.44) 

j^7^8y5y7 ^ j5j^6^4y6y2 ^ X^Z^TW^){1 - X^Z^T^^^il - X^Y^Z'^)''^ 

Notice however that in applications to tensor products, it seems we never encounter 
rational exponents but simply negative ones. The procedure is thus usually simpler. 
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7.^. Multiple su{2) products from the vector-basis construction 



A simple and different application of the formalism just developed is furnished by 
the analysis of su{2) quadruple tensor products. This application is different in that it 
does not rely on the triangle description and as such, its formulation is less direct. It 
will serve as a preparation the somewhat more complicated sp{4) example treated in the 
following section. 

The Diophantine description of this problem has been presented in section 4.2. It is 
based on the two inequalities (4.11) which are readily transformed into equalities by the 
introduction of two positive integers ai, 02: 

Ai = ni2 + ai Ai + Till - ni2 = mi2 + a2 (7-45) 

However this system does not contain any reference to the variable mn and for this reason 
we introduce the further constraint mn > which calls for a new integer variable: 

mil = as (7.46) 

We have thus a total of 8 variables : {Ai, nn, ni2, wii, TO12, ^i, ^2, ^3} and 3 equations. 
There are thus 5 independent variables, chosen to be {ai, 02, as, 7112, mi2}. The basis 
vectors, with components ordered as follows 

(ai, a2, as, ni2, mi2; Ai, nn, wii) (7.47) 

are obtained by successively setting equal to 1 one of {ai, a2, as, ni2, mi2} and the others 
equal to 0. These basis vectors together with their exponentiated version written in terms 
of appropriate grading variables read: 

(1,0, 0,0,0; 1,-1,0) :LiiV-/Ai 

(0,1, 0,0,0; 0,1,0) :N^^A2 

(0,0, 1,0,0; 0,0,1) : Mil As (7-48) 

(0,0, 0,1,0; 1,0,0) :LiiVi2 

(0,0, 0,0,1; 0,1,0) :iViiMi2 



This does not mean however that there are no diagrammatic representations for the quadruple 
product. In fact, having a set of inequahties, we can transform then into equalities, as it is 
done below, and from them set up a diagrammatic representation. In the present case, it could 
correspond to two adjacent su{2) triangles, one upside down, with their adjacent sides forced to 
be equal. Here we simply mean that the analysis will not rely on such a description. 



52 



The desired generating function is obtained from the projection to positive powers of A^n 
of the function 

(1 - LiN-'Ai){l - iVnA2)(l - LiiVi2)(l - iVnMi2)(l - M^As) ^^'^^^ 

The projection operation is done by the famihar method and the result, after setting aU 
= 1 is 

G= 1 - LiiVnMi2 

(1 - LiiVi2)(l - LiMi2)(l - Li)(l - iVnMi2)(l - iVn)(l - Mn) ^ ' ^ 

from which we read of the 6 elementary couplings Ei, ■ ■ ■ , Eq (in the order where they 
appear in the denominator) given in (4.12) and the relation E3E4 = E2E5. The above 
function is exactly the one derived in section 4.2. 

7.5. sp(4) diamonds and the vector-basis derivation of the generating function 



the integers a^, we get: @ 



.1) pertainin 


g to sp{4) 


: by setting 


ri/2 = 7^ 


Ai = 


P12 + ai 


A2 = 


TZi + a2 


^1 = 


qi2 + a5 


^2 = 


TZ2 + as 


1^1 = 


ai + ay 


1^2 = 


a4 + as 


(12 +P12 = 


as + qi2 


as + = 


a4 + T^2 


as + 27^2 = 


ag + 27^l 


ae + qi2 = 


ay +P12 



(7.51) 



This leads to a diamond-type graphical representation of the tensor product that has the 
advantage over the one presented in [18] of being linear (the sum of two diamonds is also 
a diamond): 



The original idea of looking for a diagrammatic representation of sp(4) tensor products along 
theses hnes is due to M. Walton. 
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In this picture, all data pertaining to the first (second) Dynkin label appear at the 



left (right). Dotted lines relate those two points that compose the label indicated beside 



it. Opposite continuous lines are constrained to be equal, with the length of a line being 



defined as the sum of its extremal points except for the lines delimited by the points 



(a6,7^i) and (05,7^2) where the point TZi is counted twice (the little bar besides TZi and 



TZ2 being a reminder of this particularity). Explicitly, for those lines, we have thus the 



constraint aQ+2TZi = a^-\-2TZ2. Given a triple sp(4) product, the number of such diamonds 



that can be drawn with non- negative entries yields the multiplicity of the product. For 



instance, the two diamonds that describe the triple coupling (1, 1) ® (1, 1) ® (2, 0) are: 
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The dimension of the vector basis is 8 (18 variables and 10 equations, the last four 



equations above being linearly independent). As our free variables, we choose the set 
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{TZi,TZ2,Pi2, Qi2, oi, as, og, as}- The 8 basis vectors in terms of grading variables are 

El : L2M^N2A4AlRi E2 : M-^M2N-^A~^A-^R2 



E3 : LiL^^N-''A-^A^^Pi2 E^ : L2M1N1A2A7Q12 
E5 : LiNiAi Eq : L2N2A2ASA4 



(7.52) 



E7 : MiNiA^AeAr Es : M2N2AS 

The generating function is obtained by first projecting of the function — Ei)~^ to 
positive powers for each grading variables and then by setting all grading variables equal 
to 1 except for Li, Mi, Ni^s. The sp{4) elementary couplings elementary couplings are 
simple products of the EiS (the following ^1,2,3 should not be confused with the above 
grading variables): 

Ai = E'j A2 = E^ A^ = E3E4 

Bi = Es B2 = Eq S3 = E1E2 

2 (7.53) 
Ci = E4 C2 = E2E^EqE'y C3 = EiE^E'Y 

D\ = E2E^EqE^ D2 = El D3 = E2EqE^ 

The complete list of sp(4) elementary couplings (6.2) are thus recovered. 



8. Summary and conclusion 

8.1. Listing the methods considered 

Here is the list of methods that have been reviewed here for construction tensor- 
product generating functions: 

1- The character method 

This has been extensively exemplified in section 2. It starts from first principles (the 
multiplication of two characters or more precisely, two character generating functions) but 
it is computationally complicated. 

2- Composition of generating function 

This technique uses other, simpler generating functions to construct the required 
generating function. The method has been illustrated in sections 2.3 and 2.5. One example 
of this method is the use of the Giambelli formula. 
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3- Diophantine inequalities - Elementary couplings and relations 

In this approach, we reformulate the problem of computing tensor products in terms of 
a set of Diophantine inequalities (section 4). We then look for the elementary solutions and 
their relations and from these, we construct the generating function. We can distinguish 
two ways of obtaining the elementary solutions and their relations: 

i) MacMahon's method which yields the generating function as an output of his method 
for finding the elementary solutions and their relations (here these are the syzygies of first, 
second and possibly higher orders). This approach is plagued by the technical difficulties 
of the intermediate projection operations. It is discussed in section 3. 

ii) Apply the Huet's algorithm to find the elementary couplings and then apply Grobner 
basis techniques to calculate a complete set of forbidden couplings. This is certainly 
powerful and it appears to be the most effective approach. Grobner bases are presented 
in section 5 and this method is illustrated in section 4 and 6. 

4- Diophantine equalities - The vector basis 

In this approach, we first re-express the inequalities in terms of equalities and then 
write the vector basis, relaxing the positivity constraint inherent to the elementary cou- 
plings. Since the basis vectors are independent the initial generating function is simply a 
product of terms of the form 1/(1 — E) for suitable monomials E. The elementary cou- 
plings are then recovered from the projection to positive solutions (and a second projection 
onto terms with integral exponents may also be necessary in general) and the result of the 
projection(s) is the desired generating function. This method and various examples are 
worked out in section 7. 

8.2. Conclusion 

As it has already been stressed in the introduction, the main purpose of this work is to 
prepare the ground for the analysis of fusion rules, which is the subject of a sequel paper. In 
the present work, we have reviewed the existing techniques for computing tensor-product 
generating functions and presented a comparative critical assessment of their respective 
virtues and limitations. In addition, we have focused on a model formulation linking 
generating functions to Poincare series, an idea that has first been introduced in [4] and 
further extended in [5] . Our contribution in that respect has been to rephrase this program 
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more explicitly, clarify some issues and to exemplify the procedure with many examples, 
some of which are new. 
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